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Integrable Hamiltonian systems 
on Lie groups: Kowalewski type 



By V. JURDJEVIC 



Introduction 

The contributions of Sophya Kowalewski to the integrabihty theory of the 
equations for the heavy top extend to a larger class of Hamiltonian systems 
on Lie groups; this paper explains these extensions, and along the way re- 
veals further geometric significance of her work in the theory of elliptic curves. 
Specifically, in this paper we shall be concerned with the solutions of the fol- 
lowing differential system in six variables /ii, /12, /13, Hi, H2, H3 

H2H3 f — - — ) + h2a3 - h2.a2 , 



dt \ C3 C2 

^"^^ H1H3 ( — -—)+ h^ai - hia-s 



dt \ Cl C3 

-f^i-f^2 ( 1 + hia2 - /i2ai 



dt \ C2 Cl 

dhi 



dt 
dh2 
IT 
dhs 

dt 



h2H3 


h3H2 
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C2 






Cl 


C3 


hiH2 


h2Hi 


C2 


Cl 



+ k{H2a3 - ^302) 
+ k{H3ai - FiOs) 



in which ai, 02, 03, ci, C2, C3 and k are constants. The preceding system of 
equations can also be written more compactly 

(i) = H X Q + h X a, ^ = h x Q + k{H x d) 

with X denoting the vector product in R'^ and with 

Ml 

H = I H2 \ , n = \ ^ \ , h = I h2 \ and 





When k = the preceding equations formally coincide with the equations 
of the motions of a rigid body around its fixed point in the presence of the 
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gravitational force, known as the heavy top in the hterature on mechanics. 
In this context, the constants 01,02,03 correspond to the principal moments 
of inertia of the body, while 0,1,0,2,03 correspond to the coordinates of the 
center of mass of the body relative to an orthonormal frame fixed on the body, 
known as the moving frame. The vector Q corresponds to the angular velocity 
of the body measured relative to the moving frame. That is, if R{t) denotes 
the orthogonal matrix describing the coordinates of the moving frame with 
respect to a fixed orthonormal frame, then 

^ = R{t) ^z{t) . 

The vector H corresponds to the angular momentum of the body, related 
to the angular velocity by the classic formulas = fij, i = 1,3,3. Finally, 
the vector h{t) corresponds to the movements of the vertical unit vector ob- 

served from the moving body and is given by, h{t) = R ^{t) \ . Therefore, 

solutions of equations (i) corresponding to A; = 0, and further restricted to 
/if + + /i| = 1 coincide with all possible movements of the heavy top. 

Rather than studying the foregoing differential system in R^, as is com- 
monly done in the literature of the heavy top, we shall consider it instead as 
a Hamiltonian system on the group of motions £'3 of a Euclidean space 
corresponding to the Hamiltonian function 

(ii) H = I + ^ + ^) + ai/ii + 02/12 + da/ia- 

2 V Cl 02 03 J 

This Hamiltonian system has its origins in a famous paper of Kirchhoff of 

1859 concerning the equilibrium configurations of an elastic rod, in which he 
likened the basic equations of the rod to the equations of the heavy top. His 
observation has since been known as the kinetic analogue of the elastic rod. 
According to Kirchhoff an elastic rod is modeled by a curve j{t) in a Euclidean 
space together with an orthonormal frame defined along 7(t) and adapted 
to the curve in a prescribed manner. The usual assumptions are that the rod 
is inextensible, and therefore ||^|| = 1, and that the first leg of the frame 
coincides with the tangent vector In this context, 7(t) corresponds to the 
central line of the rod, and the frame along 7 measures the amount of bending 
and twisting of the rod relative to a standard reference frame defined by the 
unstressed state of the rod. 
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Denote by R{t) the relation of the frame along 7 to the reference frame; 
then R{t) is a curve in S03(R) and therefore 



for functions ui{t),U2{t),U3{t). In the literature on elasticity these functions 
are called strains. Kirchhoff's model for the equilibrium configurations of 
the rod subject to the prescribed boundary conditions, consisting of the ter- 
minal positions of the rod and its initial and final frame, postulates that 
the equilibrium configurations minimize the total elastic energy of the rod 
^ /o^('^i^i(^) ^'^'^'ii^) + C3^'.3(t))di with ci, C2, C3 constants, determined by the 
physical characteristics of the rod, with T equal to the length of the rod. 

Prom the geometric point of view each configuration of the rod is a curve 
in the frame bundle of given by the following differential system 



with constants ai , 02 , 03 describing the relation of the tangent vector to the 
frame along 7. The preceding differential system has a natural interpretation 
as a differential system in the group of motions = k S03(M). The 
Hamiltonian H given above appears as a necessary condition of optimality for 
the variational problem of Kirchhoff. 

In contrast to the traditional view of applied mathematics influenced by 
Kirchhoff, in which the elastic problem is likened to the heavy top, we shall 
show that the analogy goes the other way; the heavy top is like the elastic prob- 
lem and much of the understanding of the integrability of its basic equations 
is gained through this analogy. To begin with, the elastic problem, depen- 
dent only on the Riemannian structure of the ambient space extends to other 
Riemannian spaces. In particular, for spaces of constant curvature, the frame 
bundle is identified with the isometry group, and the parameter k that appears 
in the above differential system coincides with their curvature. In this paper 
we shall concentrate on A; = 0, k = ±1. The case k = 1, called the elliptic 
case, corresponds to the sphere = S04(M)/S03(M), while k = —1, called 
the hyperbolic case, corresponds to the hyperboloid = S0(1, 3)/S03(M). 
As will be shown subsequently, differential systems described by (i) correspond 
to the projections of Hamiltonian differential equations on the Lie algebra of G 
generated by the Hamiltonian H in (ii), with G any of the groups E^, S04(M) 
and S0(1, 3) as the isometry groups of the above symmetric spaces. 

It may be relevant to observe that equations (iii) reduce to Serret-Frenet 
equations for a curve 7 when U2 = 0. Then ui{t) is the torsion of r(t) of 7 while 




(iii) 



lit 
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Us is its curvature K(t). Hence the elastic energy of 7 becomes a functional 
of its geometric invariants. In particular, the variational problem attached to 
J(^(K^(t) + T'^(t))dt was considered by P. Griffith a natural candidate for the 
elastic energy of a curve. Equations (i) then correspond to this variational 
problem when 02 = 03 = 1, ci = C3 = 1, and C2 = 00. 

With these physical and geometric origins in mind we shall refer to this 
class of Hamiltonian systems as elastic, and refer to the projections of the 
integral curves of the corresponding Hamiltonian vector field on the under- 
lying symmetric space as clastic curves. Returning now to our earlier claim 
that much of the geometry of the heavy top is clarified through the elastic 
problem, we note that, in contrast to the heavy top, the elastic problem is 
a left-invariant variational problem on G, and consequently always has five 
independent integrals of motion. 

These integrals of motion are H itself, two Casimir integrals 

-I- h ■ H = hiHi + /i2-f^2 + H^^h^, and two additional integrals 

due to left-invariant symmetry determined by the rank of the Lie algebra of 
G. This observation alone clarifies the integrability theory of the heavy top as 
it demonstrates that the existence of a fourth integral for differential system 
(i) is sufficient for its complete integrability. 

It turns out that completely integrable cases for the elastic problem occur 
under the same conditions as in the case of the heavy top. In particular, we 
have the following cases: 

(1) a = 0. Then, both \\h\\ and \\H\\ are integrals of motion. This case cor- 
responds to Euler's top. The elastic curves are the projections of the ex- 
tremal curves in the intersection of energy ellipsoid H = 2\~^ '''' ''f 
with the momentum sphere M = + i^l + H^- 

(2) C2 = C3 and 02 = 03 = 0. In this case Hi is also an integral of motion. This 
case corresponds to Lagrange 's top. Its equations are treated in complete 
detail in [8]. 

(3) ci = C2 = C3. Then Hiai + H2a2 + H^a^ is an integral of motion. This 
integral is also well-known in the literature of the heavy top. The corre- 
sponding equations are integrated by means of elliptic functions similar to 
the case of Lagrange, which partly accounts for its undistinguished place 
in the hierarchy of integrable tops. 

The remaining, and the most fascinating integrable case was discovered 
by S. Kowalewski in her famous paper of 1889 under the conditions that 
ci = C2 = 2c3 and 03 = 0. It turns out that the extra integral of motion 
exists under the same conditions for the elastic problem, and is equal to 

\z'^ — a{w — ka)]"^ 
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with z = 2(^1 + iH2), w = hi + ih2 and a = ai + ia2- This integral formally 
coincides with that found by Kowalewski only for fc = 0. 

The present paper is essentially devoted to this case. We shall show that 
Kowalewski's method of integration extends to the elastic problem with only 
minor modifications and leads to hyperelliptic differential equations on Abelian 
varieties on the Lie algebra of G. Faced with the "mysterious change of vari- 
ables" in Kowalewski's paper, whose mathematical nature was never properly 
explained in the literature of the heavy top, we discovered simple and direct 
proofs of the main steps that not only clarify Kowalewski's method but also 
identify Hamiltonian systems as an important ingredient of the theory of el- 
liptic functions. 

As a byproduct this paper offers an elementary proof of Euler's results of 
1765 concerning the solutions of 

dx dy 

/PR /P(y) 

with P an arbitrary fourth degree polynomial with complex coefficients. Com- 
bined further with A. Weil's interpretations of Euler's results in terms of addi- 
tion formulas for curves = P{x), these results form a theoretic base required 
for the integration of the extremal equations. 

This seemingly unexpected connection between Kowalewski, Euler and 
Weil is easily explained as follows: 

The elastic problem generates a polynomial P{x) of degree four and two 
forms R{x, y) and R{x, y) each of degree four satisfying the following relations 

(iv) R{x,x) = P{x), and R\x,y) + {x - yfR{x,y) = P{x)P{y). 

We begin our investigations with these relations associated to an arbitrary 
polynomial P{x) = A + 4Bx + 6Cx^ + ADx^ + Ex^. In particular, we explicitly 
calculate the coefficients of R corresponding to R{x,y) = A + 2B{x + y) + 
3C(a;^ -f- y"^) -\- 2Dxy{x + y) + Ex^y^. Having obtained the expression for R, 
we have Re{x,y) = —{x — y)^9'^ + 2R(x,y)9 + R{x,y) is the form in (iv) that 
corresponds to the most general form Re{x,y) = R{x,y) — 9{x — y)^, that 
satisfies R{x,x) = P{x). 

We then show that Ro(x,y) = contains all solutions of f^, — = 

as 9 varies over all complex numbers. This demonstration recovers the result of 
Euler in 1765, and also identifies the parametrizing variable 9 with the points 
on the canonical cubic elliptic curve 

with 



92 



= AE -4BD + 3C'^ and 53 



= ACE + 2BCD- AD'^ - B'^E -C^ 
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via the relation 9 = 2(^ + C). The constants g2 and are known as the 
covariant invariants of the elhptic curve C = {{x,y) : = P{x)}. Andre Weil 
points out in [13] that the results of Euler have algebraic interpretations that 
may be used to define an algebraic group structure on F U C The "mysteri- 
ous" change of variables in the paper of Kowalcwski is nothing more than the 
transformation from C xC into T x F given by {N, M) e C x C to (Oi, O2) in 
F X F with Oi = iV - M and O2 = + M. 

The actual formulas that appear in the paper of Kowalewski 
d^i dx dy ^ d^2 dx dy 

are the infinitesimal versions of Weil's addition formulas. 

Oddly enough, Kowalewski omits any explanation concerning the origins 
and the use of the above formulas, although it seems very likely that the 
connections with the work of Euler were known to her at that time (possibly 
through her association with Weierstrass) . 

The organization of this paper is as follows: Section I contains a self- 
contained treatment of Hamiltonian systems on Lie groups. This material 
provides a theoretic base for differential equations (i) and their conservation 
laws. In contrast to the traditional treatment of this subject matter geared to 
the applications in mechanics, the present treatment emphasizes the geometric 
nature of the subject seen through the left-invariant realization of the sym- 
plectic form on T*G, the latter considered as G x g* via the left-translations. 

Section II contains the reductions in differential equations through the 
conservation laws (integrals of motion) leading to the fundamental relations 
that appear in the paper of Kowalewski. Section III contains the proof of 
Euler's result along with its algebraic interpretations by A. Weil. Section IV 
explains the procedure for integrating the differential equations by quadrature 
leading up to the famous hyperelliptic curve of Kowalewski. 

Section V is devoted to complex extensions of differential system (i) in 
which the time variable is also considered complex. Motivated by the brilliant 
observation of Kowalewski that integrable cases of the heavy top are integrated 
by means of elliptic and hyperelliptic integrals and that, therefore, the solutions 
are meromorphic functions of complex time, we investigate the cases of elastic 
equations that admit purely meromorphic solutions on at least an open subset 
of under the assumption that ci = C2, while the remaining coefficient C3 
is arbitrary. We confirm Kowalewski's claim in this more general setting that 
the only cases that admit such meromorphic solutions are the ones already 
described in our introduction. In doing so we are unfortunately obliged to 
make an additional assumption (that is likely inessential) concerning the order 
of poles in solutions. This assumption is necessitated by a gap in Kowalewski's 
original paper, first noticed by A.A. Markov, that apparently still remains open 
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in the literature on the heavy top. We conclude the paper with an integrable 
(in the sense of Liouville) elastic case that falls outside of the meromorphically 
integrable class suggesting further limitations of Kowalewski's methods in the 
classification of completely integrable elastic systems. 



1. Hamiltonian systems on Lie groups 

We shall use g to denote the Lie algebra of a Lie group G, while g* will 
denote the dual of g. The cotangent bundle T*G will be identified with G x g* 
via the left-translations: an element {g,p) in G x g* is identified with ^ G T*G 
by p = dL*g^ with dL* denoting the pull-back of the left-translation Lg{x) = gx. 

The tangent bundle of T*G is identified with TGxg* xg*, the latter further 
identified with G x g* x g x g*. Relative to this decomposition, vector fields on 
T*G will be denoted by (X[g,p),Y*{g,p)) with {g,p) denoting the base point 
in T*G and X and Y* denoting their values in g and g* respectively. Then 
the canonical symplectic form u on T*G in the aforementioned trivialization 
of r*G is given by: 

(1) a;(,,p) ((Xi,y;), {X2,Y,*)) = Y*{X,) - Y*{X2) - p[X^,X2] . 

The correct signs in this expression depend on the particular choice of the Lie 
bracket. For the above choice of signs, [X, = Y{Xf) — X{Yf) for any 

function /. 

The symplectic form sets up a correspondence between functions H on 
T*G and vector fields H given by 

(2) u^g^p){H{g,p),V) = dH^g^p){V) 

for all tangent vectors V at {g,p). It is customary to call H a Hamiltonian 

function, or simply a Hamiltonian, and H the Hamiltonian vector field of H. A 
Hamiltonian H is called left-invariant if it is invariant under left-translations, 
which is equivalent to saying that H is a function on g*; that is, H is constant 
over the fiber above each point p in g*. 

For each left-invariant Hamiltonian H, dH, being a linear function over 
g*, is an element of g at each point p in g*. Then it follows from (1) and (2) 
that the Hamiltonian vector field if of a left-invariant Hamiltonian is given by 
{Xip),Y*ip) ingxg* with 

(3) X{p) = dHp and Y*{p) = -ad* {dHp){p). 

In this expression ad*X denotes the dual mapping of adX : g — given by 
adX{Y) = [X,Y]. 
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It immediately follows from (3) that the integral curves {g{t),p{t)) of H 
satisfy 

(4) dLg-i(t)^ = dHp^t) and ^ = -ad*(di?p(t))(p(i)) , 
and consequently 

(5) Pit) = Ad*(,)(p(0)) 

with Ad* equal to the co-adjoint action of G on q*. Thus, the projections 
of integral curves of left-invariant Hamiltonian vector fields evolve on the co- 
adjoint orbits of G. 

When the group G is semisimple the Killing form is nondegenerate and 
can be used to identify elements in g* with elements in q. This correspondence 
identifies each curve p{t) in g* with a curve U{t) in g. For integral curves of a 
left-invariant Hamiltonian H, the equation ^ = —ad* {dHp){p{t)) corresponds 
to 

(6) ^ = 

The expression (6) is often called the Lax-pair form in the literature on the 
Hamiltonian systems. 

We shall use {F, H} to denote the Poisson bracket of functions F and H. 
Recall that {F, H}{g,p) = uj(^g p'j{F{g), H(p)}. It follows immediately from (1) 
that for left-invariant Hamiltonians F and H, their Poisson bracket is given 
by {F, H}{p) = p{[dFp, dHp]), for all p in g*. 

A function F on. T*G is called an integral of motion for H ii F is constant 
along each integral curve of H, or equivalently if {F, H} = 0. A given Hamil- 
tonian is said to be completely integrable if there exist n — 1 independent in- 
tegrals of motion Fi,... , Fn-i that together with F^ = H satisfy {Fi, Fj} = 
for all The independence of Fi,... is taken in the sense that the 
differentials dFi, . . . ,dFn are independent at all points on T*G. 

Any vector field X on G lifts to a function Fx on T*G defined by Fx{0 = 
^{X{g)) for any ^ G T*G. In the left-invariant representation G x g*, left 
invariant vector fields lift to linear functions on g* , while right-invariant vector 
fields lift to Fx{g,p) = p{dLg-i o dRgX^) with denoting the value of X at 
the group identity e of G. The preceding expression for Fx can also be written 
as Fx{g,p) = Ad*-i{p){Xg). Therefore, along each integral curve {g{t),p{t)) 

of a left-invariant Hamiltonian H, 

Fx{g{t),p{t)) = Ad*_i(,)p(t)Xe 

= Ad*.i(,)oAd;(,)(p(0))Xe = (p(0))Xe 

and consequently Fx is an integral of motion for H. 
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The maximum number of right-invariant vector fields that pairwise com- 
mute with each other is equal to the rank of g. Consequently, a left-invariant 
Hamiltonian H always has r-independent integrals of motion Poisson commut- 
ing with each other, and of course commuting with H, with r equal to the rank 
of 0. 

In addition to these integrals of motion, there may be functions on q* that 
are constant on co-adjoint orbits of G. Such functions are called Casimir func- 
tions, and they are integrals of motion for any left-invariant Hamiltonian H. 
On semisimple Lie groups Casimir functions always exist as can be seen from 
the Lax-pair representation (6). They are the coefficients of the characteristic 
polynomial of U{t) (realized as a curve on the space of matrices via the adjoint 
representation) . 

With these concepts and this notation at our disposal we shall take g 
to be any six dimensional Lie algebra with a basis Bi, B2, B3, Ai, A2, A3 that 
satisfies the following Lie bracket table: 
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The reader may easily verify that the following six dimensional matrices 
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satisfy the above Lie bracket table under the matrix commutator bracket 
[M, N] = NM - MN. For A; = 0, is the semi-direct product ik sosiR), for 
= 1, is so4(i?,), and for k = —1, g = so(l, 3). 

Throughout this paper we shall use hi and Hi, to denote the linear func- 
tions on Q* given by hi{p) = p{Bi), and Hi{p) = p{Ai), i = 1,2,3. These 
functions are Hamiltonian lifts of left-invariant vector fields induced by the 
above basis in q. Finally, as stated earlier, we shall consider a fixed Hamilton- 
ian function H on q* given by 

1 fHf Hi Hl\ 

2 V ci C2 C3 y 

for some constants ci, 02,03 and 01,02,03. 

We shall refer to the integral curves of H as the extremal curves. For each 
extremal curve {g{t),p{t)), x{t) = g{t)ei will be called an elastic curve. Elastic 
curves are the projections of the extremal curves on the underlying symmetric 
space G/K with K denoting the group that stabilizes ei in (written as 
the column vector, with the action coinciding with the matrix multiplication). 
It can be easily verified that K ~ ^0^{R) and that G/K is equal to R^,5^ 
or depending whether k = 0, 1 or —1. The remaining columns of g give 
the coordinates of the moving frame vi{t),V2{t),v^{t) defined along x{t), and 
adapted to the curve x{t) so that = a\Vi{t) + a2V2{t) + 03^3 (i). 

The semi-simple case. For /c / 0, the Killing form T is nondegenerate 
and invariant in the sense that T{[A, B],C) = T{A,[B,C]). We shall take 
T{A,B) = itrace(A5). It follows that T{A,B) = -(E-^iOjOi + kbibi), with 
A = Yh^i ^i^i + ^i^i ai^d ^ = Yh=i "'iA + hBi. Upon identifying p in g* 
with U in Q via the trace form, we get that 
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and equation (6) yields the following differential system: 
(8) 



dhi 


h2H3 


h3H2 ^ 
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"dt 




C2 
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"dt 


Cl 
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H1H3 
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Cl 


C3 
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_ H1H2 


H1H2 


dt 


C2 


Cl 



+ (/i2a3 - ^302) , 
+ {hia2 - h2ai) . 



Remark. The foregoing differential equation can also be obtained by using 
the Poisson bracket through the formulas 

^ = {hi,H}, and = {Hi,H}, i = l,2,3. 

Apart from the vector product representation given by (i) of the introduc- 
tion, differential system (8) has several other representations. The most imme- 
diate, that will be useful for Section V, is the representation in so3(i?) x so3(i?) 

obtained by identifying vectors ^4 = 0:2 I in with antisymmetric matri- 

\«3/ 

03 02 \ 

—ai I . In this representation differential system (8) 
0/ 

dP 

[n,K] + [A,P], — = [n,P] + k[A,K] 



^Hi \ /hi\ I ai 

^^^2 , P = U2 , and 1 = as 

The characteristic polynomial of the matrix U in (7) is given by 

\'' + \W\Hf + k\\hf) + {H-hf ■ 

hence, 

(10) K2 = \\h\\'^ + k\\H\\'^ and K3 = h ■ H 
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are the Casimir functions on g. Being constant on each co-adjoint orbit of 
G, they Poisson commute with any function on q*, and in particular they 
Poisson commute with each other. Since g is of rank 2, it follows that in 
addition to K2 and K3 there are two extra integrals of motion for H by our 
preceding observations about right-invariant vector fields. Together with H 
these functions constitute five independent integrals of motion, all Poisson 
commuting with each other. So H will be completely integrable just in case 
when there is one more independent integral that Poisson commutes with H. 

The Euclidean case. The group of motions is not semisimple, hence 

the Hamiltonian equations cannot be written in the Lax-pair form as in (6). 
The following bilinear (but not invariant) form reveals the connections with 
the equations for the heavy top: 



{A,B) = y^^ajQi + bjbj, 



1=1 



with 



A = ^ aiAi + hiBi and -B = ^ aiAi + hiBi. 



i=l 



Relative to this form every p = Y17=i ^i-^i + ^i^t 5* identified with 
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/o 










hi 







I] 









-Hi 
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-H2 


Hi 
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m fl. 



Then along an extremal curve {g{t),U{t)) of H, functions FL{g,U) 

1 

X R 

terms of the translation x and the rotation R, we see that Fl becomes a func- 
tion of the variables x, R, h, H and is given by 



{U,g ^Lg) are constant for each L in g. Upon expressing g 



m 



FL{x,R,h,H) = h-{R-^{v + Vx) + H -R-^V 



with 



L 
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Vl 







V2 









-Vl 




-V2 


Vl 


/ 



and V 




The Lie algebra of E^ is of rank 3, because all translations commute. Taking 
L in the space of translations amounts to taking y = 0, and so functions 



F^{x,R,h,H) = h-R- 
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Poisson commute with each other, and are also integrals of motion for any 
left-invariant Hamiltonian. The functions Fy form a three dimensional space 
with Fi = Rh • ei, F2 = Rh -62, F^, = Rh ■ 63 a basis for such a space. The 
elements of this basis are not functionally independent because of the following 
relation: 

f2 + f| + f| = hi + hi + hi 

Consequently, the functions Fy give at most two independent integrals of mo- 
tion. 

By using the Poisson bracket. Table 1 one shows that the differential 
equations for U{t) can be written as 

= h{t) X h{t), ^ = H{t) X h{t) + h{t) X a. 

Hence, the Hamiltonian equations in this case coincide with equations (8) 
and (9) for k = 0. Consequently, the conservation laws defined by (10) apply 
to = and we get that 

K2 = \\h{t)\\'^ = constant, and = H(t)-h{t) = constant 

along the integral curves of H. 

Together with H, K2 , -f^s and any two functions among Fi , F2 , F3 account 
for five independent integrals of motion all in involution with each other, and 
the question of complete integrability in this case also reduces to finding one 
extra integral of motion. 

Having shown that the equations (i) in the introduction coincide with 
the Hamiltonian equations (9) we now turn to the integrable cases. Since the 
extra integrals of motion for the three cases mentioned in the introduction are 
evident, we shall go directly to the case discovered by Kowalewski. 

So assume that ci = C2 = 2c3 and that 03 = 0. Normalize the constants 
so that c = C2 = 2 and C3 = 1. Then, equations (8) become: 

^ = H^h2 - ^-^2^3 - ka2H2., = ^-ff2i^3 - 02/13 , 

^ = \Hihz - Hzhi + kaiHs, = —HiHs + mhs, 

at 2 dt 2 

^ = ^iH2hi - Hih2) + k{a2Hi - aiH2), = 02/^1 - ai/i2- 

Set z = ^{Hi + iH2), w = hi + ih2, and a = ai + ia2- Then 
— = —{H^z-ahs) and — = i[h3Z - Hsw + kHsa). 
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Let q = z — a{w — ka). Then, 

dq dz dw -itt i ^ \ ■ n tt , tt \ 
— = 2z— (^~T' — "^(-"3-2^ ~ an^z) — la^h^z — li^w + kH^a) 



= -iHi{z^ - aw + ka^) = -iH2.{t)q{t) . 



Denoting by q the complex conjugate of q we get that ^ = iHs{t)q{t), and 
hence 

—q{t)q{t) = -iH2,qq + iH2,qq = 0. 

Thus, 

q{t)q{t) = \q{t)^ = constant. 

Hence — a{w — ak)\'^ is the required integral of motion. 
We shall refer to this case as the Kowalewski case. 



2. The Kowalewski case: Reductions and eliminations 

It will be convenient to rescale the coordinates so that the constant a is 
reduced to 1. Let 

a f t \ l^^/t\ a f t \ l,/t\ 

^3 



^ = ^ rr ' ^3 = 77-^3 (ti ] ' y = vri^ rr ' ^3 = 77^ 

It follows from the previous page that 

(11) ^ = -^{x3X-y3), ^ = iivsx - xsy + kxs) , 

— = Imy, and = (Im + 2A;Im x) 

dt dt 

are the extremal equations in our new coordinates. 

The integrals of motion in these coordinates become: 

(12) H = }-{Hl + Hl) + ^Hl + aihi + a2h2 

= + -i^l + Reaw = |a|^(xa; + -x| + Rey) , 

K2 = \\hf + k\\Hf = \a\'^{\y\'^ + yl + k{4\x\^ +xl)) , 
K3 = h-H = |ap(2Rexy + xsys) , 
Kl = qq = \aWx^-iy-k)f 

and therefore we may assume that a = 1. This rescaling reveals that system 
(11) is invariant under the involution 



'^{x,y,X3,y3) = {x,y,-X3,-y3). 
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We shall now assume that the constants if, K2,Kz, K4 are fixed, and use V 
to denote the manifold defined by equations (12). Now F is a two dimensional 
real variety, contained in R^, that can be conveniently parametrized by one 
complex variable according to the following theorem. 

Theorem 1. V is contained in the set of all complex numbers x and q 
that satisfy 

(13) P{x)q + P{x)q + Ri{x,x)+Kj{x-xf = 

with 

P{x) = k2-2K2.x + 2Hx^ - x'^, and 

(14) 

Ri{x,x) = {HK2 - Kl) + 2K^k{x + x) + {2Hk - 3K2){x^ + x^) 
+ 2K3Xx{x + x)- Hx^x"^ + {Hk - 2K2){x - xf , 

where H = 2H - 2k, k2 = K2 - kH - K^. 

Proof. Equation (13) is a consequence of eliminating X3 and from the 

constraints (12) as follows: 

Begin by expressing the integrals of motion in terms of x,x,q,q, X3 and 
ys- Putting y = x'^ — q + k, in the expression for H leads to 

H = xx+ ^xl + Re(x^ -q + k) = xx + ^xl + ^(x^ + x'^ - {q + q) + 2k). 

This relation simphfies to 2H — 2k = H = {x + x)^ — {q + q) + x^. 

Then, K2 = x'^x^ + A;^ + j/| + kH + 2kxx — {x'^q + xP'q) and hence, 

K2-kH -k'^ = K2 = x^x^ + yi + 2kxx - {x^q + x'^q). 

Finally, = (xx + k)(x + x) — {xq + xq) + xzyz- 

Eliminating y^ and X3 from the preceding relations leads to: 

(14a) 

{K^ - {xx + k){x + x) + {xq + xq)f' = xlyl 

= {K2 + x'^q + x^q - x^x^ - 2kxx){H - {q + q)- {x + xf). 

The homogeneous terms of degree two in q and q in the preceding expres- 
sion reduce to 

{xq + xq)"^ — {xP'q + x^q){q + q) 

2-2 I -2 2 I o 2-2 - /-2 2 , 2 - , -2 - , 2-2\ 

= X q + x q + 2x x qq — {x q +xqq + xqq + xq) 
= -Ki{x-x)^. 

Therefore, relation (14a) can be reduced to 

Pq + Pq + R = 
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for suitable polynomials R and P in the variables x and x. It follows that 

P = {K2- x^x? - 2kxx) + x'^{H - {x + xf) - 2{K3 - {x + x){xx + k))x 
= K2- 2K3X + x^{H - 2k) -x^- x^x^ - 2kxx - x^x^ - 2x^x 

+ 2kxx + 2x^x^ + 2x^x 
= k2-2Kzx + 2Hx^ -x'^. 

Thus P is a polynomial of degree 4 in the variable x only. 
Then, x) = Ri{x, x) + Kl{x — x)^ with 

= {H-{x + xf){K2-x'^x'^-2kxx)-{K3-{x + x){xx + k)f. 

The expression for Ri further simplifies to 

Ri = {Hk2- Kl) + 2kK3{x + x) + {2Hk-3K2){x'^ + x'^) + 2K3Xx{x + x) 
- H^x'^x'^ + {Hk - 2K2){x - xf 

by a straightforward calculation. 

Equation (13) identifies x as the pivotal variable, in terms of which the 
extremal equations can be integrated by quadrature. For then q is the solution 
of (13), and the remaining variables are given by 

X3 = H — [x + x) + {q + q), and 7/3 = K2 — x x + x q + x q — 2kxx. 

Theorem 2. Each extremal curve x{t) satisfies the following differential 
equation: 

(15) = p(.x)+Qm^-^?, 

with P{x) as in the previous theorem. 

Proof. It follows from equation (11) that —4 (^)^ = {x^x — y^f. 

{X3X - y^f = xlx^ - 2x3y3X + y| = {H + (q + q) - (x + x)'^)x^ 

— 2x(i^3 — (x + x)(x.T + /c) + ixq + xq)) 
+ {K2 — 2kxx — x^xP' + .x^g + xP'q) 

= K2- 2K3X + ^x^ + (x^(g + ^) - (x + x)^x^ 
+ 2x(x + x)(xx + A;) 

— 2x(xg + XQ') — 2A;xx — x^x^ + x^q + sP'q) . 

But then, 

x^{q + q) — {x + x)^x^ + 2x(x + x)(xx + k) — 2x{xq + xq) 
— 2A;xx — x^x^ + x^q + x^g = q{x — x)^ — x^ + 2/ex^, 
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and therefore 

(xsx - y^f = K2- 2K3X + Hx'^ + 2kx^ - x^ + q{x - xf 
= P{x) + q{x — xf'. 

Theorem 3. Let Ro{x,x) = K2 - Ks{x + x) + H{x^ + x^) — x'^x^ — 
{H-k){x-xf. Then 

(16) Ro{x,x) = P{x), and B^{x,x) + {x - xfRi{x,x) = P{x)P{x) 
where R\ has the same meaning as in Theorem 1. 

Proof. Let C, = x^x — y^. We shall first show that Rq = 

CC = {X3X - y3){x3X - ys) = x^xx - X3ys{x + x) + y^ 

= {H +{q + q)- {x + xf)xx + (^2 - x'^x^ - 2kxx + x'^q + x^q) 

— {K3 + xq + xq — {x + x){xx + k)){x + x) 

= K2 — i^3(x + x) + k{x + xf — 2kxx + Hxx — x^x^ 
+ {xx{q + q) — xx{x + xf + x'^q + x^q 

— {xq + xq){x + x) + xx{x + xf). 

The above expression reduces to 

K2 - Ksix + x) + k{x'^ + x^) + Hxx - x^x^ 

because 

+ ^) + x^g + x^q — {xq + xq){x + x) = xx{q + q) — xx{q + = 0. 

Since 2k + H = 2H, the preceding expression can also be written as 

K2 - Ks{x + x) + H{x'^ + y^) - x^x^ - {H - k){x - xf 

showing that Rq = CC- 

It follows from the proof of Theorem 1 that = P{x) + q{x — xf. There- 
fore, i?o(x,x) = P(x) and, 

Rl{x,x) = iCCf = C'C' = iP{x) + q{x - xf){P{x) + q{x - xf) 

= P{x)P{x) + P{x)q{x - x)^ + P{x)q{x - x)^ + qq{x - x)^. 

Thus 

RI{x,x) = P{x)P{x) - Ri{x,x){x - x)^, 

because P{x)q + P{x)q = —Ri{x,x) — K2{x — x)'^ as can be seen from relations 
(13) in Theorem 1. Our theorem is proved. 
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The essential relations obtained by the preceding calculations, can now be 
summarized for further reference as follows. 

For each choice of constants of motion, the algebraic variety defined by 
cqTiations (12) is parametrized by a single complex variable x through the 
relations 

P{x)q^ + {Ri{x,x) + Kl{x-xf)q + KlP{x) = and qq = Kl, 

with P{x) a polynomial of degree four, and Ri{x,x) the form defined by (14). 
The form i?i has a companion form Rq{x, y), also of degree four, that satisfies 

Ro{x,x) = P{x) and R^{x,x) + {x — x)^Ri{x,x) = P{x)P{x). 

Finally, —4 (^)^ = P{x) = q{x — x)'^ is the extremal differential equation 
that needs to be solved. 

Apart from the more general nature of the constants that occur in the 
foregoing expressions, all of the above relations are the same as in the original 
paper of Kowalewski from 1889. It will be convenient to refer to the above 
relations as Kowalewski's relations. 



3. Lemniscatic integrals and addition formulas of 
L. Euler and A. Weil 

We shall now show that Kowalewski's relations provide remarkable in- 
sights into the theory of elliptic curves and elliptic integrals starting with the 
very beginning of the subject with the work of G.S. Fagnano in 1718 concern- 
ing the arc of a lemniscate, and the subsequent extensions of Euler concerning 
the solutions of a more general differential equation 

where P{x) is an arbitrary fourth degree polynomial. ([12]) 

To begin with, note that when the constant of motion is equal to zero, 

then 




Thus our extremal equations already contain the cases studied by Euler in con- 
nection with the length of an arc of a lemniscate. It is therefore not surprising 
that our study must be closely related to the work of Euler. However, it is 
remarkable that the fundamental relation 

Rl{x,x) + {x-xfRi{x,x) = P{x)P{x) 



INTEGRABLE HAMILTONIAN SYSTEMS ON LIE GROUPS 



623 



of Kowalewski provides an easy access to Euler's famous results concerning the 
solutions of 

dx _^ dy ^ Q 

/PR ^/P{y) 

To explain Euler's results and its relevance to this paper we shall assume 

that 

P{x) = A + ABx + GCx'^ + ADx^ + Ex'^ 

is a general polynomial of degree four. We shall also denote by R{x,y) = 
A + 2B{x + y) + 3C(x2 + y^) + 2Dxy{x + y) + Ex'^y^ a particular form that 
satisfies R{x,x) = P{x). 

Theorem 4. 
R{x,y) = - AB'^ + 4{AD - 3BC){x + y) 

+ 2{AE + 2BD - 9C^){x'^ + y^) + 4{BE - 3CD)xy{x + y) 
- AD'^x^y'^ - {AE + ABD - 9C'^){x - yf 
is the unique form that satisfies 

R\x,y) + {x-yfR{x,y) = P{x)P{y) . 
In particular, R{x, x) defines a polynomial Q given by 

Q{x) = A{-B^ + 2{AD - 3BC)x^ + {AE + 2BD - 9C^)x'^ 
+ 2{BE -SCD)x^ - D^x"^). 

Proof. Consider F{x,y) = P{x)P{y) — R^{x,y). Because P{x) = R{x,x), 
P'{x) = 2§{x,x). Therefore, ^{x,y) = P'{x)P{y) - 2R{x,y)^{x,y), and 
for X = y 

dF 1 

— (x,x) = P'{x)Pix)-2P{x)-P'{x) = 0. 

It follows by an analogous argument that ^{x,y) = ioT y = x. Thus 
both F and its partial derivations ^ and ^ vanish at a; = y. Consequently, 

F{x,y) = (x — y)'^R{x,y) for some binary form R{x,y). 

Now, by differentiation of y) + {x — y)'^R{x, y) = P{x)P{y), 



dR OR 
2R{x, y)^{x, y) + 2(x - y)R{x,y) + {x - yf — {x,y) = P\x)P{y), 



and 



fdRdR d^R\ ^ , ,dR , ,dR 



oyox 
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Let Q{x) denote R{x, x), and set a; = y in the above expression. It follows 

that 

which upon substituting 2^ (x, x) = P'(x) reduces to 

P'{xf + 2P{x)^{x,x) - 2Q{x) = P\xf. 
I oyax 

Solving for Q gives 

^^^^ =^^^)S^^'^^"i^'^^^' 

= {A + ABx + QCx^ + ADx^ + Ex^){mx + AEx^) 

- ^{4B + 12Cx + 12Dx'^ + AEx^f 

= -4:B^ + {8AD - 2ABC)x + (32BD + AAE - 2ABD - 36C^)x^ 

+ {ASCD + 125^ - 72CD - 8BE)x^ 

+ {32D^ + 2ACE - 361)2 _ 24EC)x'^ 
= -AB^ + 2{AAD - 12BC)x + {8BD + AAE - 36C^)x'^ 

+ 2{ABE -\2CD)x^ - AD'^x^ 
= A' + AB'x + eC'x^ + AD'x^ + E'x"^. 
Therefore, 

R{x, y) = A' + 2B'{x + y) + -t{x^ + y^) + 25xy + 2D'xy{x + y) + E'x'^y^ 

with 27 + 25 = 6C", where 7 and 5 are to be determined from the fundamental 
relation B? + {x — y)^R = P{x)P{y). Upon equating the homogeneous terms 
of degree 4 in the fundamental relation we get: 

9C2(a;2 + y^f + 8BDxy{x + yf + 2AEx^y'^ - 7(x - yf{x^ + y^) + 25xy{x - yf 

= AEy"^ + IQBDxy^ + mC'^x^y'^ + 16BDx^y + AEx"^. 
It follows that 

AE{x + yf -9C'^{x + yf + 8BDxy = -/{x'^ + y'^) + 25xy 
and therefore 

7 = AE-gC"^, and 6 = AE + ABD -90"^. 
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The substitution of these values in the expression for R leads to 

R(x, y) = - + {AAD - l2BC){x + y) + {ABD + 2AE - lSC'^){x^ + y^) 

+ {ABE - l2CD)xy{x + y) - AD^x'^y'^ 

- {AE + ABD-9C^){x-yf 
and the proof of the theorem is finished. 

Let Ro{x,y) = R{x,y) — 9{x — y)^ denote the most general biquadratic 
form that satisfies Rg{x,x) = P{x) with 9 an arbitrary parameter. Then 
y) = — (x — + 2R{x, y)9 + R{x, y) is the unique form that satisfies 

(17) Rl{x,y) + {x-yf^e{x,y) = P{x)P{y). 

Remark. For the relations of Kowalewski obtained earlier in the paper, 
P = K2 — 2K'iX + 2Hx^ — x^, and the forms Rq and Ri are given by Rq = Rq, 
with ^ = - /c, and i?i = -{x - yf{H - kf + 2R{x, y){H - k) + R. 

Now, to return to the general case, $g is symmetric with respect to x and 
y and can be written as 

^e{x, y) = ae{x)y'^ + 2he{x)y + ce{x) = ae{y)x'^ + 2be{y)x + co{y) 

for some quadratic expressions a0,b0 and cq. 

Writing R = a + 2/3(x + y) + S'y{x'^ + y^) + 2Sxy{x + y) + ex^y'^ - C,{x - yf 
we have 

^0 = -{x- y)\9^ + C) + 2{A + 2B{x + y) + ^C{x^ + y^) 
+ 2Dxy{x + y) + Ex^y^ + {R + ({x - yf), 

and therefore, 

ae{x) = {2E9 + e)x^ + {AD9 + 25)x + {<oC9 + 87 - (0^ + C)) , 
ce{x) = {6C9 + 37 - {9"^ + C))^^ + {4B9 + 2/3)x + {2A9 + a) , 
beix) = {2D9 + 6)x^ + {9^ + C)x + {2B9 + (3) . 

After substitution of the values for a, l3, 7, 6, e, ( given by Theorem 4, the pre- 
ceding expressions become 

ae{x) = {2E9 - AD'^)x^ + {AD9 + A{BE - 3CD)x + AE - {9 - 3Cf , 
C0{x) = {AE -{9- 3Cf)x^ + {AB9 + A{AD - 3BC))x + {2A9 - AB^) , 
bg{x) = {2D9 + 2{BE - 3CD))x^ + (^^ - QC^ + AE + ABD)x 
+ {2B9 + 2AD -3BC). 

Theorem 5. Let Gq denote the discriminant bg — aecg. Then Gq{x) = 

p{9)P{x) with 

(18) p{9) = 29{9 - 3Cf + 29{ABD - AE) + AB'^E + AAD^ - 2ABCD. 
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Proof. Let X be any root of P{x). Then, 
{R{x, y) - e{x - yff + {x - yf{-{x - yfO^ + 2R{x, y)e + R{x, y)) = 0. 
When 

^e{x,y) = 0, R{x,y) - 9{x - yf = 

and therefore ^e{x, y) = has a double root. Therefore = 0. This argument 
shows that P{x) is a factor of Gq. Since Gg is a polynomial of degree four in 
x, it follows that Ge{x) = p{6)P{x) for some polynomial p{9)- It now follows 
by an easy calculation, using the explicit expressions of ag, bg and cq, that p{d) 
is given by expression (18). This ends the proof. 

Expression (18) is most naturally linked with the cubic elliptic curve T = 
{(^) v) '■ ff = — (72^ — 53}, where 52 and g2 are the invariants of C = {(x, u) : 
rt^ = P{x)} exphcitly given as follows: 

(19) 52 = AS-4SL> + 3C^ = ACE + 2BCD - AD"^ - B'^E - C^. 

The identification is obtained as follows: Let = 2{^ + C). Then, p{^) = 
4(4^3 _ _ 4^£, ^ 3(^2)^ _ ^ 2BCD - AD^ - B'^E - C^)). 

Letting ry^ = | we obtain rj^ = 4^^ — ~ 53- 

The next theorem is a paraphrase of the classical results of Euler ([3]). 

Theorem 6. ^g(x,y) = is a solution for —^S= ± fy = for each 

number 9. Conversely, for every solution y{x) of either differential equation 
f"^, lb f^, = there exists a number 9 such that ^gix^y) = is equal to 

Proof. Consider ^g{x,y) = for an arbitrary number 9. Then 

^ ^ -bojy) + p^GM y ^ -hejx) +aVG^) 

agiy) ag{x) 

with Gg = b"^ — agCff, and p = ±1, a = ±1. Hence 

1 9^0 

(21) 2~dx' = ^"^(y) + ^^*^2/) = PvG'e(y) and 

1 9^0 , 

2 "ay" = yae(a:^) +?'e(a;) = o-VGe(a::) • 

Using the results of Theorem 5 we get that the curve yix) of (20) satisfies: 



dy ^ Id^ ^ -p^Gsiy) ^ -pyJp{9)P{y) 

dx dx I dy a^/G^ a^p{9P{x) ' 

and therefore 

dy p dx 

+ - , = 0. 
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Thus, v(x) is a solution of — H — ^= = when ^ = 1. Otherwise, it 
is a solution of ^ p= = . 

Suppose now that yfx) is a particular solution of either — ;= — h = 

or i^— M= = 0. Let x = a, y = bhe any point such that 6 = via). We 

need to show that there exists 9 such that $6t(a, 6) = 0, or that 



Since 

we get that 
(22) 



R{a, b) ± yR'^ia, 6) + (a - b^Ria, b) 
^ " (a - 6)2 ■ 

ii:^(a,6) + (a-6)^^(a,6) = P{a)P{b) 



R{a, b) ± ^P{a)P(b) 



(6 -a) 



2 



The appropriate sign for ^ depends on whether y{x) is a solution of J^^— + 

= 0, or ^= ?==• The correct way of choosing the sign for 9 

\/P{y) v^(^) v^(^) 

will be made clear through the discussion of the related addition formulas of 
A. Weil. 

Addition formulas of A. Weil. A. Weil points out in his paper [13] that 
the results of Euler (Theorem 6) have algebraic interpretations that may be 
used to define a group structure on F U C Recall that 

r = {(e,r?): = 4^' - 52^ - 53} and C = {{x,u) = = P{x)}. 

It turns out, quite remarkably, that these algebraic observations of Weil 

provide easy explanations for some of the formulas used by Kowalewski, and 
for that reason it will be necessary to explain Weil's interpretation of Euler's 
results. 

With each solution ^g{x,y) = Weil associates two transformations 
from C into C, depending whether they change ^ into or into — each 
parametrized by the points of T rather than 6. More precisely, let P = (^, rj) 
be any point of F. Let ^0 denote the form corresponding to ^ = 2(,^ + C). For 
each point M = {x,u) of C define two points iVi = (yi,fi) and N2 = {y2,V2) 
on C by the following formulas: 

—b0(x) + 2r]u 1 

(23(a)) yi = -— , vi = - — (xae(yi) + 60(2/1)) , 

ag{x) 2ri 

(23(b)) y2 = -^^(^) -^^" , y, = _J_(^a,(y2) +6,(2/2)) • 
ae{x) 2r] 
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It follows that 

dyi _ 8^0 d^0 _ 2r]vi dy2 _ d^e ,9^6 _ -"^^2 

dx dx dy 2r)u ' dx dx dy 2r)u ' 

consequently, the mapping (P, M) — Ni takes ^ into ^ , while the mapping 
(P, M) — > 7^2 takes ^ into 

Following Weil we shall write P + M = A'"i and P- M = N2, with the 
understanding that — P = — r/) and — Af = (,t,— ?7,). Wc assume that on T 
the group law coincides with the usual group law, with the point at infinity on 
r acting as the group identity. 

Formula (22) may be used to show that for each M and N on C there 
exist points P and P' on F such that 

P + M = N and P' - M = N 

by the following simple argument: let P = (^,r/) and P' = Note that 

P' = 2M when N = M, while P is the point at infinity when N = M. Thus, 
^ = 2(6 + C) and ^' = 2(9' + C) with 9 and 9' appropriately chosen among 

R{x, y) ± uv 
(y - x)2 

For uv = sj P{x)P{y), '^^^'^^-j gives a finite value when x = y because 

R{x,y)-uv ^ R{x,y)-^P{x)P{y) 
{x - yf {x - y)2 

^ (P(x,y) - ^P{x)P{y)){R{x y) + ^P{x)P{y)) 

[x - yy{R{x, y) + ,JP{x)P{y)) 
^ R^{x,y)-P{x)P{y) ^ ~R{x y) 

[x - yf {R{x, y) + ^P{x)P{y)) R{x, y) + ^P{x)P{y) ' 

and the latter expression evidently tends to ~2p{^'^ when y tends to x. Thus, 

^, ^ R(x, y) - uv 
(x - yY 

Then use formula (23(b)) to define 77' = —^{xae>{y) + bo'{y)). 

The other choice of sign uv = — y^P{x)P{y) leads to 9' = by an 

analogous argument. In such a case define rj' = ^{xaQi{y) + bo{y)). 

The values of r] corresponding to 9 are given by 

ri = -^{xa0{y) + 60 (y)) when uv = ^/P{x)P{y), and 
■q = l-{xa0{y) + 50 (y)) when uv = ~^JP{x)P{y). 
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To these addition formulas of Weil we can add their infinitesimal versions 
that appeared in the paper of Kowalewski. 

Theorem 7. Let O = rj) and O' = (^', -q') denote points of T, and 
let M = (x^u) and N = {y,v) denote points of C related by the formulas 
= N-M, 0' = N + M. Then, 

d$' dx dy , d£ dx dy 

= — + — , and — = + — . 

rf u V rj u V 

Proof Since ^e{x,y) = -{x - y)'^6^ + 2R{x,y)9 + R{x,y), 

dMx,y) = {-2{x-yfe + 2R)de + —^dx + —^dy. 

ox oy 

For {x, y, 9) for which ^e{x, y) = 0, d^g{x, y) = 0, and 9 = ^^(^^^J^i" depending 
whether uv = -s/ P{x)P{y), or uv = —y^P{x)P{y). 

Assuming that uv = \J P{x)P{y), we have 9 = ^(^'^)^y7" and = 
xag{y) + bg{y) = -2ryf and = yae{x)+bg{x) = 2r]u. Thus, 

—2uvd9 — Arjvdx + Arjudy = 0, 

and since d9 = 2d^, 

d$, dx dy 

r] u V 



In the remaining case uv = —y/P{x)P{y), and = 2riv, and = 2-qu. 

Thus, again 

d^ dx dy 

rj u V 

The corresponding differential form for ^' is obtained by analogous argu- 
ments and its derivation will be omitted. 

With these formulas behind us we finally come to the integration of the 
extremal differential equations. 

4. Kowalewski's integration procedure 

In her original paper of 1889 Kowalewski made the following change of 
variables: 

Ro{xr,X2)- ^P{XI)P{X2) , 1„ 
= 2{x,-X2f + 2^^ 

and 



R0{XI,X2) + ^P{XI)P{X2) , 1„ 

2{xi-X2Y 2 
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and then claimed that 

dsi dxi dx2 ^ ds2 dxi dx2 

where S[s) = 4s^ — g2S — g^. 

In her notation, i?o(a;i, X2) = —XiX2 + 6iiXiX2 + 2ico{xi+X2) + CQ — k'^, 
and P{x) = -a;^ + 6£ix'^ + 4£cox + eg - A:^. Then, 

i?o(a;i,X2) = -x\x\ + Mi{x\ + X2) + 2£co{xi + 2:2) + Cq - /c^ - 3£i{xi - ^2)^ 
= i?(a;i,X2) - 3^i(xi - 0:2)^ 

and so Kowalewski's constants are related to the constants of this paper as 
follows: 

Sii = H, cl-k^ = K2, 2ico = -K3 . 

The change of variables, which has often been called "mysterious" in the 
subsequent literature, is nothing more than the mapping from C xC into F x F 
that assigns to M = {xi,ui) and N = {x2,U2) the values Pi = = 
M — N and P2 = {s2, S{s2)) = M + N . The variables si and S2 in her paper 
coincide with ^1 and ^2 in this paper by the following calculation: 



R{xi,X2) - Mi{xl - xl) - ^P{xi)P{X2) , 1, 
2[xi-X2y 2 



_ R{xi,X2)-^yP{x)P{y) 3^1 1 
2(xi-X2)2 2 2 

_0i . _ 2(^1 + C) _ 

-y-^i - — ^ ii - 6 

because £1 = C. An analogous calculation shows that S2 = ^2- Therefore, 
Kowalewski's change of variables is fundamentally a paraphrase of Theorem 7 
in this paper. 

The extremal equations in terms of the new coordinates s\ and S2 are 
obtained exactly in the same manner as in the original paper of Kowalewski. 
For convenience to the reader we include these calculations. 

For the sake of continuity with the rest of the paper we shall use 
and (y, f ) to denote the points of C, rather than and (x2,U2), and use 

{ii,rii) and {S,2,'>]2) instead of {si,S{si)) and {s2, S{s2))- Upon identifying y 
with X, we get that 

= P{y) + m{x{t)-y{t)f. 
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Since ^ = (— l)*i7 + ^ with i = 1,2, it follows that along each extremal 
curve x(t) 

y u2 1^ ^ + ^2 V y ^ ' uv dt dt ■ 



The substitution of 



A\dt 



2 



^= + and 

Ro{x{t),y{t)) (Theorem 2) 
into the preceding expression leads to 

_ 2^2^^ + (x - yf{v^q + u^q) + 2(-l)^i?o(x, y)™ 



2u^v^ - (x - + - yY) + 2(-l)*i?o(a;, 



^2^2 



because v^g + it^g + i?i + — y)^ = (Theorem 1). 

With the aid of R\{x — y)^ = — iilg + u^v^ the above expression simplifies 

to 

4 (diiV v?v'^ + Rl + 2{-iyRouv-{x-yYKl 



Tjf \ dt J u^v^ 

_ {uv + i-lfRoY-ix-yYKi 



A{x-yY (Ro + {-lfuvV Kl 



y2y2 I 2(x — yY ) 4 y 

Recall that ^2{x~-y)^^ ~ ^« ~ + 1 ■ Upon introducing two new constants 
H k Ki , , H k Ki 



into the above equality we get 

1 [diiV {x-yY 



rjf \ dt J u^v^ 
Finally, ^2 - 6 = 1^^^ and hence 
diY\^ _ U{^i) ^^^^^ ^^^^^ ^ _^^^3 _ ^^^^ _ ^^^^^^ _ ^^^^^^ _ 
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Consequently, 



+ p — ;^^= = U With p = ±1. 



The last equation coincides with the celebrated equation of Kowalewski 
when p = 1. 



5. The ratio of coefficients and meromorphic solutions 

Kowalewski begins her investigations of the heavy top with a brilliant 
observation that both the top of Euler and the top of Lagrange, being integrable 
by means of elliptic integrals, admit Laurent series solutions of complex time 
for their general solution. In the papers listed in the bibliography she classifies 
all the cases of the heavy top whose solutions are meromorphic solutions of 
complex time, concluding that meromorphic solutions require that two of the 
coefficients must be equal to each other, say ci = C2, in which case 

( 1°) Cl = C2 = C3, and 01,02,03 arbitrary, 
( 2°) Cl = C2, C3 arbitrary, but oi = 02 = 0, 
( 3°) Cl = C2 = 2c3, and 03 = 

are the only situations admitting general meromorphic solutions. 

In this section we shall pursue Kowalcwski's approach in a more general 
setting and classify all the cases in which the Hamiltonian solutions of 

1 fHf Hi Hl\ ^ 

i/ = - — + — + — + aihi + 02/12 + 03/13 

2 V Cl C2 C3 y 

are meromorphic functions of time under the assumption that ci = C2. 

We shall also include the following limiting cases in the subsequent analy- 
sis: 

H = -Hi + oi/ii + 02/12 + 03/13 

obtained as the common value ci = C2 tends to 00 while C3 = 1, and 

H = ^{hI + Hi) + aihi + a2h2 + azhz 

obtained as C3 tends to 00 while ci = C2 = 1. 

These limiting Hamiltonians are important for different reasons. The first 
case is completely integrable but falls outside of Kowalewski's classification, 
while the second case, which also falls outside of Kowalewski's classification, 
arises naturally in the theory of curves. 
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In what follows, we shall assume that all the variables Hi, H2, H^,hi,h2, /13 
are complex (but the coefficients ci, 02,03 and 01,02,03 are real) and search 
for the cases where the solutions of 

(24) ^ = [^^K] + [a,Pi ^ = [n,P]+k[a,K] 

are meromorphic functions of complex time t, for at least an open set of initial 
values in C^. 

It will be convenient to use c instead of C3, and to use m to denote the 
common ratio ^ = ^. That is, m = :r- and hence ci = ^ = C2- Therefore, 

Cl C2 ' Cl ^ m ' 

n = I ™i?2 I . It follows that JqKJq = O with 





The limiting case H = + aihi + 02/12 + 03^3 occurs when m = and 

c = 1, while the other limiting case H = \{H\ + i?|) + a\hi + 02/12 + 03/13, 
occurs when ^ = 1 and lime = 00, 

We shall use {A,B) to denote the trace of —\AB for any antisymmetric 
matrices A, B with complex entries (not to be confused with the Hermitian 
product on C^). Then, 

{A,B) = A B = 0161 + 0262 + 0363 ■ 
In terms of this notation 

H = ^{JoKJo,K) + {a,P). 
It is easy to verify that the Casimir functions 

G = (P, P) + k{K, K) and J = (P, K) 

remain constants of motion for the complex system (24). 

We shall now seek conditions on the ratio m and the coefficients oi, 02, 03 
so that solutions K{t),P{t) of equation (24) are of the form 

K{t) = r^'^iKo + Kit + • • • + K„t" + •••), 

P{t) =t-^^PQ + p,t + --- + p^e + ---) 

for some positive integers ni and n2, and matrices Kn,Pn in so3(C) for each 
n = 0, 1, . . . , with neither Kq nor Pq equal to zero. 

Evidently not every solution has a pole in C. For instance, K = and 
P = a is a solution for any choice of the constants in question. We shall 
assume that the meromorphic solutions occur for at least some open set of 
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initial conditions in C , which in turn imphes that the meromorphic solutions 
should be parametrized by six arbitrary complex constants. 

Lemma 1. ni = 1 and n2 = 2, provided that [JqKqJq, Pq] ^ 0. 

Proof. The leading terms in ^ = [fi,P] + k[a,K] are given by 
-n2t~("2+i)pQ = t-("2+"i)[JoKoJo,Po] which implies that n2 + 1 = ni + n2. 
Therefore, ni = 1. 

Then the leading terms in ^ = [n, K] + [a, P] are given by 

-t-^Ko = t~^[JoKoJo,Ko]+t-''^[a,Po]. 

If [a,Po] = then —Kq = [JoKqJo, Kq]. The latter relation can hold only for 
Kq = 0. Thus [a,Po] 7^ and therefore n2 = 2. The proof is now finished. 

In the original paper Kowalewski claims the results of Lemma 1 without 
any proof ("on s'assurc facilement, en comparant les exposants des premiers 
termes dans les membres gauches et dans les membres droits des equations 
considerees que Ton doit avoir ni = I, mi = 2," [9, p. 178]). This claim, 
first criticized by A. A. Markov, is still considered an open gap in the origi- 
nal approach of Kowalewski. The assumption [JqKqJq, Pq] 7^ is very likely 
inessential (for the contrary would reduce the number of possible solutions), 
and hence the original investigations of Kowalewski are in all probability cor- 
rect in spite of the gap. We shall subsequently assume that n\ = l and n2 = 2 
without going into details caused by vanishing of the above Lie brackets. 

Upon equating the coefficients that correspond to the same powers of t in 
equations (24) we come to the following relations: 

(25) -Ko = [JoKoJo,Ko] + [a,Po], 2Po = [JoKoJq, Pq] 

and 

n 

(26a) (n - l)Kn = ^[Joi^i Jo, Kn-i] + [a,Pn] , 

1=0 
n 

(26b) (n - 2)Pn = J2[JoKiJo, Pn-i] + k[a, Kn-2] 

1=0 

for n>l. 

The same procedure applied to the Hamiltonian H and the Casimir func- 
tions G and J gives 

1 " 

(27) -^{JoKiJo, Kn-i) + {a,Pn) = d2nH , 

1=0 



(28) 



n n—2 

J2{Pi,Pn-i)+kJ2{Ki,Kn-2-i) = S4nG, 
1=0 1=0 
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(29) ^{Pi,Kn-i) = 53nJ , 



1=0 



with Sij denoting the Dirac function equal to 1 only for i = j, and otherwise 
equal to zero. It follows that 



Y,{JoKiJo,K2-i) + {a,P2) = H 

i=0 

-1 2 

Y,{Pi,P4-i) + kY,{Ki,K2-i) = G 

i=0 1=0 
3 



1=0 

and consequently the first four stages of equations (26) must generate three 
arbitrary constants H, G and J. 

Following Kowalewski's original paper we shall also assume, without any 
loss in generality, that 02 = and that a\ + a\ / 0. This situation can 
always be arranged by a suitable rotation of a. We shall also rule out the cases 
ai = as it corresponds to the case of Lagrange, and also rule out m = 1 as it 
corresponds to the well-known case ci = C2 = C3. 

We shall now use PniQni^n to denote the entries of while fn^9nihn 
will denote the entries of P„ in the recursive relations (25) and (26). We then 
have: 

Theorem 1. There are finitely many solutions of equation (25) if and 
only if 2m — 1 7^ 0. They are given by 
(a) 

2eic 2c 
Po = 0, qo = , ro = 0, fo = isho, go = 0, ho 



m m(az + eiai) 



and 



(b) Po = -i^go, qoi'^rn - 1) = vq = 2eic, 

ai 

fo = — , 90 = isfo, ho = 
ai 

with = 1 in either (a) or (b). 

When 2m — 1 = 0, then 03 = and go is an arbitrary complex number 
corresponding to each case e = ±1. 
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Proof. The equation 2Po = [Po, JqKoJo] is linear in Pq and can be written 

as 

/ 2 -Iro -fqo\ (fo\ 
(i) hro 2 J ^ J = . 

V?^o -?;>o 2 y V/^o/ 

The determinant of the preceding matrix is 2(4+ ^(pg + Q'o) + ^''o)' ^"^^ 
must vanish to admit nonzero solutions. Hence, 



(ii) '^{pl + qo) + ^ro = -4 . 



It also follows from 2Po = [Pq, JoKqJo] that {Po,Po) = and that {Pq, JqKqJq) 
= 0. Therefore, 

(iii) /o + 9o + h-l = 0, and m(po/o + qo9o) + ro/io = . 

Casimir relation (29) implies that {Pq, Kq) = Po/o + 9050 + ^0^^0 = 0, which 
together with (iii) implies that r^h^ = 0. Consider first the case tq = 0. 
The relation Kq = [Kq., JqKqJq] + [Pq, a] gives 

(fYl — \\ [TTl — 1) 

(iv) Pq = roqo-asgQ, qq = ropo + oa/o -oi/io, = oi^o • 

c c 

Hence 50 = 0, and consequently /g + = 0. Thus, /o = eihQ with £ = ±1. 
It follows from (27) that \^{pl + ql) + + ai/o + 03/10 = 0. Thus, 

/ . N, r , I'TT'/ 2 2n -"i^c(pn + 9n) 4c 2c 
(aiez + a3)/io = ai/o + as/^o = -^^(Po + Qo) = ^2 = ^ = ^ 

(m 7^ 0, as can be seen from (ii)). Therefore, 

2c 



ho 



m(a3 + eioa) 
Equations (iv) now imply that pq = 0, and 

qo = asfo - aihQ = a^eho - oi/iq = ho{-ai + £^3) 

2 , , 2eic 

-(— ai + £za3) 



m(a3 + eiai) m 

We have now verified the solutions given by (a) . The solutions in (b) correspond 
to /iQ = 0, and are obtained as follows: 

Linear system (i) implies that ^(/oQc — goPo) = 0- Combined with (hi) 
this relation gives that either m = 0, or that p^ + Qq = 0. In both cases 

= —4 as can be seen from (ii). Therefore, vq = 2eic with e = ±1. 

Then go = = ^ (equations (iv)), and /o = = ff, because 

^Tq + ai/o + 03/10 = 0. That is, gQ = eifo. 
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The relation pofo + goffo = gives fo{po + ^igo) = 0. Thus po = -eiqo- 
Then, 

go = J'oPo + as/o = (2e^)c(-£^go) + 03/0, 

c c 

or 

go(l + 2(m-l)) = aa/o = 
Hence solutions in (b) are also verified concluding the proof of the theorem. 

It remains to resolve the recursive relations (26). They are linear in 

and Pn and can be expressed as 
(30) 

n—1 

(n - l)Kn + [Kn, JqKoJq] + [Ko, JoK^Jq] + [Pn, a] = ^[Joi^i Jo, K^-i-i] , 

n—1 

(n - 2)P„ + [P„, Joi^o Jo] + [Po, Joi^n Jo] = X^t-^oi^i Jo, J'n-l-i] + k[a, i^„_i] 

i=l 

forn > 1 with = 0. 

This linear system in six variables PmQ.mfmfmgni hn can be written more 
explicitly as follows: 

(n - l)pn + ^{mgnro - go^n + rnqoVn - qnVo) - a^Qn = An , 

{n - l)qn + ^{rnPo - mpnTo + roPn - mpoVn) + azfn - aiK = Bn , 
(n - l)rn + aign = Cn 

and 

(n - 2)fn + - {mqnho - goVn + mqoK - gnVo) = , 
c 

(n - 2)gn + -(rn/o - mpnho + ro/„ - mpohn) = , 

Tfi 

{n -2)hn-\ {Pngo - foqn + PoSn " /n^o) = Fn 

c 

with An, Bn,Cn, Dn, En, Fn denoting the appropriate quantities on the right- 
hand side of equation (30). 

1. Meromorphic solutions for tq = 0. It follows from Lemma 2 that 
Po = ro = go = 0. Therefore the preceding equations cascade into the following 
independent subsystems: 

(31) (n - l)pn H qoVn - as^n = A , 

(n - l)rn + aiQn = Cn , 

{n-2)gn + -{rnfo-mpnho) = En 



638 



V. JURDJEVIC 



and 

(32) (n - l)qn + a^fn - aiK = , 

{n-2)fn + -{mqohn + mqnho) = , 

Tfli 

(n - 2)hn - — {hqn + Qofn) = Fn ■ 

The determinant A of the overall system is the product of determinants 
Ai and A2 with 



Then, 



/ 1 (m— 1) 

/ n - 1 ^—^qo 



Ai = det 







n- 1 



-03 \ 
ai 



K-T^o Ifo n-2j 
/ n - 1 03 -ai \ 



A9 = det 



^ho n-2 ^qo 



m 



m , 



m 



A2 = (n - l)[(n - 2)2 + —q^] ho{{n - 2)a^ - ai— go) 



m 



m 



— /o(a3— go + ai(n - 2)) 



m 



m 



(n - l)((n - 2)2 + — g2) _ (n - 2)-(a3^o + /ooi) 



TTl 

+ -^go(ai/io - 03/0). 

RecalUng that go = /o = «£^o and ho = ^i^as+ieai) S^* 
A2 = (n- l)((n-2)2 -4) -2(n-2) +4 = (n + l)(n - 2)(n - 4). 
Then, 

Ai = n - I n - 1 (n - 2 - ^/q - -/io ^ -—qo + asin - 1 

c c c c 

= (n - l)((n - l)(n - 2) - ai— - —^003) ^ ai/iogo 

= (n — l)((n — l)(n — 2) ) -{m{m — l)goai + (n — l)(m — 1)03) 



= (n- l)((n- l)(n-2) 



''0 



m 



(m — l)(2£iai — (n — 1)03) . 
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But 
ho 



{2eiai — {n — 1)03) = 



2{2eiai — (n — 1)03) 
m(a3 + eiai) 



-(2ai — (n — l)a| + ^0103(1 + n)). 



For Ai to be equal to zero the imaginary part of Ai must be equal to 
zero, which occurs only when 03 = 0. But then 

Ai ={n- l)((n - l)(n - 2) - -) - ^^"^ ~ 

m m 

_ (n - l)((n - l)(n - 2)m - 2) - 4(m - 1) 
m 

_ (n - 3)((n2 - n + 2)m - 2) 

m 

Therefore, A = when n = 2, n = 4, and n = 3 provided that 03 = 0. 
In such a case A also vanishes for another positive integer value k given by 
{k^-k + 2)m-2 = 0. 

Theorem 2. All solutions of equations (26) are parametrized by at most 
four arbitrary constants, and hence do not provide for general solutions of 
equation (24). 

Proof. The maximum number of solutions of equations (26) occurs when 
n = 3 is a zero of A, that is when 03 = 0. Our theorem follows from the simple 
observation that the kernel of the cascaded linear system with 03 = 0, 



(n - l)p„ 



m 



, (n - l)qn - aihr. 



, 



{n - l)r„ + aign = , 



1 Tfl/ 

{n - 2)fn + -mqohn H qnho = , 

c c 



1 TTl TTl 

{n-2)gn + -{rnfo-rahopn) = 0, (n - 2)/i„ /og„ go/n = 

c c c 

is at most one dimensional at each singular value of A. Recall that 

2c 2eic 
fo = eiho, ho = and that go = ■ 

mai m 

The foregoing linear systems are described by the following matrices: 
/ n — 1 — ai \ /n — 1 2eim{m — 1) \ 



— n 

ai 



^ mai 



- 2 2£i 
2ei n-2) 



and 







V -- 

^ ai 



n 



2ei 
01 



ai 
n-2/ 
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The rank of each of these matrices is at least two, and therefore there 
can be at most four constants arising from n = 2, n = 3, n = 4 and 
(A;^ — k + 2)m — 2 = 0. Our proof is finished because the number of parameters 
required for general solutions of equation (24) is six. 



Remark. It might seem plausible that equations (26) do not admit solu- 
tions at the singular stages. Remarkably, that is not the case. In fact, for 
n = 2 





with 



/2 = f^?2 + ^fli) ^''^d /i2 = with q2 an arbitrary complex number. Then, 
for n = 3 

/ ei(m— 1) \ 




, with /i3 = --kaiq2 



\ -f93 J 
and an arbitrary complex number. Finally, for n = 4 



K4 




Pa 



2^ ^2 - 2ei«4 




with ^4 arbitrary. 

2. Meromorphic solutions for ho = 0. In this case it is convenient to take 
c = 1. Then the solutions of equations (25) are given by 

po = -eiqo, qo{2m - 1) = — , tq = lei, fo = — , go = i£fo, ho = 



with e = ±1. We shall first suppose that 2m — 1 ^ and show that the number 
of constants that parametrize the general solution is less than six. 

Since po + isqo = and fo + sigo = it is natural to consider the following 
change of variables: 



— Pn "I" ^^Qm "^n — Pn ^^Qm "^n — fn "I" ^^Qm — fn isqn- 
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(n — l)un — ie{m — l)rou„ — ie{m — l)uorn — iea^Wn — aiiehn 

= Un{n - 1 + 2{m — 1)) — ieasWn — ieaihn = An + ieBn , 
Vniin ~ 1) ~ 2(m — 1)) + ie{m — l)vorn + iea^Zn + ieaihn = An — ieBn , 

/ ix , I -ix eiaiWn , eiqzn „ 
(n - l)r„ + aiQn = [n - l)r„ \ — = Cn , 

(n - 2 - 2)wn = Dn + ieEn , 

(n - 2 + 2)zn - eizQVn + eimvohn = Dn + ieEn , 

(n - 2)hn + m{pngo - foqn + Po9n - fnQo) 



= {n-2)hn + mfoeiun - qoWn 



Fn ■ 



The matrix corresponding to this linear system is given below: 



M = 



2m + n — 3 






V m/( 



n + 1 — 2m ie{m — l)vo 



o£^ 










n — 1 


-eizQ 




-iea^ 
ieai 



2 2 



n-4 



—leai \ 
ieai 







n eimvQ 
-go n-2 ./ 



The determinant A of the preceding matrix is given by A = (n — 4) (n + 
1 — 2m) Ai with Ai equal to the determinant of 



n-l ^ 



ieai \ 




/n-3+2m 


—eizo n eimvg 

\ mfoei n-2/ 
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Then 



Ai = (ra-3 + 2m)(n-2)(n(n- 1) 



- mfo£i{-ieai{n{n - 1) — )) 



aizo 



n(n - 1) - -^)((n - 3 + 2m)(n - 2) - maifo) 



n{n - 1) - 2)((n - 3 + 2m)(n - 2) - 2m) 
n(n - 1) - 2)((n - l){n - 2) + 2(m - l)(n 
n + l)(n - 2)((n - l)(n - 2) + 2(m - l)(n 
n + l)(n - 2)(n - 3)(n - 2 + 2m). 



2) 
2) 



2m) 
2m) 



Therefore, 
A 



(n + l)(n - 2)(n - 3)(n 



Theorem 3. 2m -1 = and 03 = 
meromorphic solution of equation (24). 



- 4)(n + 1 - 2m)(n - 2 + 2m). 
is the only case admitting a general 



Proof. When m 7^ and 2m — 1 7^ then there are four constants that 
parametrize the solutions corresponding to n = 2, n = 3, n = 4 and one of 
the factors n — 2 + 2m = or n + 1 — 2m = because the kernel of M is one 
dimensional for each singular vahic of A. When m, = 0, the kernel of M is two 
dimensional for n = 2 and therefore contributes two constants that together 
with the constants corresponding to n = 3 and n = 4 again account for four 
constants. 

In the remaining case 2m — 1 = 0, and the zero-th stage (Theorem 1) 
accounts for two arbitrary constants, provided that 03 = 0, that together 
with the constants corresponding to n = 1, n = 2, n = 3, n = 4 produce 
six arbitrary constants, which is the exact number required for the general 
solution. 

This analysis shows that 



(i) 


a - 


= 0, 




(ii) 


Cl 


= C2 


ai = 02 = 0, 


iii) 


Cl 


= C2 


= C3, a arbitrary. 


iv) 


Cl 


= C2 


= 2c3, 03 = 0, 



are the only cases that admit meromorphic solutions thus extending the con- 
clusions of Kowalewski to a wider class of systems. 

In the literature on Hamiltonian systems this claim of Kowalewski is often 
confused with the classification of completely integrable systems, that is, with 
systems that admit six independent integrals of motion all in involution with 
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each other. This issue was considered by R. Liouville in 1896 ([11]) where 
he claimed that there is an extra integral of motion whenever the ratio ^ 
is an even integer. The exact nature of this claim is not altogether clear as 
the arguments in Liouvillc's paper are difficult to follow. In addition, the 
homogeneity properties of the equations for the heavy top, upon which any of 
Liouville's arguments are based, do not hold in our setting when k ^ and 
further limit the relevance of Liouville's paper. However, the following simple 
example shows that there arc additional cases of completely integrable systems 
which are outside of meromorphic class and therefore cannot be integrated by 
means of Abelian integrals. 

3. A completely integrable system not of Kowalewski type. This case cor- 
responds to the limiting ratio m = and 03 = 0. Then H = -|-ai/ii + 02/12, 
and 

—7—= H2Hz-Ha2, —7— = --H1-H3 + fiaai, — j— = "'201 - /i2ai , 
dt dt dt 

^= h2H3-kH3a2, ^ = -hiHs + kHsai, = k{Hia2 - H2ai) . 

dt dt dt 

Let w = hi + ih2- Then, ^ = —iH3{w — k{ai + ia2)) ■ Hence, ^{w — ka) = 
—iHs{w — ka), with a = ai + ia2, and consequently \w — ka\'^ is a constant of 
motion. As in all other integrable cases, this case is also completely integrable 
in the sense that there are six integrals of motion all in involution with each 
other. 

Let be an angle defined hy w — ka = Re^^ . Along each solution curve, 
- ka) = Rie'^^ = -iHaiRe'^) and hence ^ = -H3. Then 

1 / d6 \ ^ 

= ^ = H- Redw = H- Red{ka + Re'^) 
= H — A;|ap — {Rai cos 6 — Ra2 sin 6) . 



2 \dt 



So even though 9 can be expressed in terms of elliptic integrals, it follows 
from the preceding results that the remaining equations cannot be integrated 
in terms of rational functions of 6. 

The other limiting case H = ^{Hf + + ai/ii + 02/12 + 03/13 may be 
symplectically transformed into H = ^{Hf + + hi, provided that 02 = 
03 = and ai = 1. H corresponds to the total elastic energy of a curve 7 
in M given by ^ J {'^'^{t) + T'^{t))dt with nit) and r(t) denoting the curvature 
and the torsion of 7. While it is not known whether this Hamiltonian system 
admits an extra integral of motion, it nevertheless follows from the foregoing 
analysis that its equations cannot be integrated on Abelian varieties in terms 
of meromorphic functions. 
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